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Abstract 
O, 

CN , It has been conjectured that the two-dimensional M = 2 Wess-Zumino model with a quasi- 

^ I homogeneous superpotential provides the Landau-Ginzburg description of the M=2 super confer mal 

minimal models. For the cubic superpotential W = A<l>'^/3, it is expected that the Wess-Zumino 
model describes A2 model and the chiral superfield ^ shows the confermal weight (h, h) = (1/6, 1/6) 
at the IR fixed point. We study this conjecture by a lattice simulation, extracting the weight from 
the finite volume scaling of the susceptibility of the scalar component in <1>. We adopt a lattice 



< 

-(— > 

^ . model with the overlap fermion, which possesses a Nicolai map and a discrete R-symmetry. We 

set a A = 0.3 and generate the scalar field configurations by solving the Nicolai map on L x L 

(N 

K^ ' lattices in the range L = 18 - 32. To solve the map, we use the Newton-Raphson algorithm with 

C^^ ■ various initial configurations. The result is 1 — h — h = 0.660 it 0.011, which is consistent with the 

''^ , conjecture within the statistical error, while a systematic error is estimated as less than 0.5 %. 
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Introduction — Two-dimensional Af = 2 Wess-Zumino (WZ) model with a quasi- 
homogeneous superpotential is expected to describe, at the infrared (IR) fixed point, the 
M = 2 superconformal minimal models iHS]. This conjecture, which followed an an alog ous 
discussion for A/" = by Zamolodochikov[7], has been tested in various aspects |8l-ll5|. If 
one could calculate correlation functions of the WZ model directly in the IR region, it would 
give us a further test of the conjecture. Since the coupling of the WZ model becomes strong 
in the IR region, this certainly requires non-perturbative techniques. 

Lattice methods can be useful for this purpose. Although the realization of supersym- 
metry on the lattice is known to be difficult because of the lack of translation invariance 

nn 

and the failure of the Leibniz rule[16|, 117|, one can preserve some part of the extended su- 



)ersymmetries which are not directly related to the translation symmetry |18|-|23| 24j | 25 1 
261]. In such a lattice model, in general, extra fine tunnings are required to keep the model 
within the universality class of the target continuum theory. But, in some lower-dimensional 
models, the part of the supersymmetries preserved on the lattice, together with other lattice 
symmetries, turn out to be sufficient to suppress the extra relevant and marginal operators. 
The two-dimensional M = 2 WZ model is an example of such supersymmetric field theories 



that can be formulated successfully on the lattice in the above sense 18|-|20| , and can be 
studied non-perturbatively through lattice simulations. 

The purpose of this article is to provide a non-perturbative numerical evidence for the 

m$) = A<I)73 by sim- 



27 



and a discrete R- 



above conjecture in the case of the simplest cubic superpotentia _ 

ulating the lattice WZ model which possesses the Nicolai map 

symmetry 20|] . For the cubic superpotential, it is expected that the WZ model describes the 

A2 model (c = 1) and the chiral superfield $ shows the conformal weight (/i, h) = (1/6, 1/6) 

at the IR fixed point. We will extract the conformal weight of the chiral superfield $ from 

the finite volume scaling of the susceptibility of the scalar component in $. Interestingly, 



the A2 model is also realized by the Gaussian model (c = 1' 
the M = 2 supersymmetric points, K = l/127r or 3 /An 3 



30, 



with the coupling constant K at 



3l| . We will extract the coupling 



constant K of the Gaussian model by identifying the phase factor of the scalar component 
(j) as the 27r-periodic Gaussian field. This will provide a clear numerical evidence of the full 
recovery of A/" = 2 super symmetry in the IR limit. 

Lattice formulation of WZ model — We adopt the lattice WZ model which possesses a 



Nicolai map and a discrete R-syinmetry 20[|. The model is formulated with overlap Dirac 



operatorldy, D = (1/a) [l+X/yJX^X] , where X = l-(l/2) [(l+7^)aV+ + (l-7M)«V;] 



with the first forward, backward difference operators Vjj. D satisfies the Ginsparg- Wilson 



relation 



34 



M' 



351], -D73 + 73i5 = with 73 = 73(1 — O-D), and it can be expressed as 
D = -yiSi + 72S'2 + aT in a spinor decomposition, defining three difference operators Si, 
S2 and T in relation of 2T = — (Si)^ — (5*2)^ + (aT)"^. Then, the bosonic and fermionic 
actions are given by 

5b = a^Y,{<P*i'^T)(j) + W*'{l-a^T/2)W' (1) 

+W'{-Si + iS2)(p + W*\-Si - 1S2W], 
Sf = a'J2^{D + F{<p))^, (2) 

X 

where F(0) = l±iiw"^ + l^W*"^ and W = dW{(f))/d(j), W" = d^W {(f)) / d(j)\ In 
the finite lattice of the volume L x L, we adopt the periodic boundary condition for both 
bosonic and fermionic fields. 

Throughout this article we consider the cubic superpotential Vr($) = A$'^/3. The cou- 
pling A, which has the mass dimension one, is the unique mass parameter of our model, 
besides the lattice spacing a. Thus A gives the scale under which the WZ model reduces to 
a conformal field theory. To see the conformal behavior on the lattice, we should prepare a 
large lattice size L/a ^ (aA)~^, while the continuum limit is L/a — )■ 00 and aX — )■ 0. 

We summarize the symmetries of our lattice model. First we note that the WZ model 

n 

possesses the Nicolai map [27|, which is exphcitly given by 

rj = W' + a{<p- ^W')T + (0* - ^W*'){Si + iS^). (3) 

The Jacobian of this map from {</>,(/)*} to {?7,?7*} precisely cancels the overlap fermion 
determinant \D + -F(0)|, while the bosonic action S-q is identical to the Gaussian weight 
a^ X]x I^(^)p5 which allows us to interpret r/, rj* as the random white noises. Thanks to this 
map, the lattice model has the following on-shell nilpotent supersymmetry Q: 

(4) 
(5) 
(6) 



(^<p = -^_, 


(^<p = -tp. 


W+ = -v\ 


Q^- = -V 


Q^+ = 0, 


Q^- = 0, 



where we write ■?/' = (■?/'i,-?/'2)*, ip = {i>i,i>2), 'ip± = {ipi±ip2) / V^ ^.nd iIj± = {ipiT'4^2)/V^- This 
lattice model also has a Z3 R-symmetry. S-p has the U{1) R-symmetry under — )■ e~^*"0, 
ip — )■ e^^'^'^ip, ip — )■ ■?/'e*"'^^. But this symmetry is broken to Z^ {a = nn/Sjii G Z) by the last 
two terms in Sb, which would be surface terms in the continuum limit. 

With these symmetries, we can show in the lattice perturbation theory that the desired 
continuum limit is achieved without extra fine-tunings. The redefined fields (/? = A0, x — -^"^ 
are helpful for this discussion. In this notation, ip has mass dimension 1 and x has 3/2, and A 
is factorized as the overall factor 1/A^ in the action. This overall factor counts the number of 
loops as the Planck constant does. Consider the generic radiative correction to an operator 
O of mass dimension p (> 0) in the action, 

''''''' ' + cia^-^ + C2d> + ■■■] [d^xO (7) 



where Cq, ci, ... are constants. The first, second and third terms represent the contributions at 
tree, one-loop and two-loop levels. In the continuum limit a — ?■ , the corrections terminate 
at the two-loop level. At the tree level the lattice action agrees with that of the WZ model 
in the continuum limit. So we have to consider the operators with p < 2. Such operators 
which preserve the Z3 R-symmetry and the fermion number are a constant and ^p*^p. But 
the constant has no effect on the path-integral and (/9*(/9 is forbidden by the supersymmetry 
Q. Thus we do not need any extra fine-tunings to achieve the desired continuum limit, at 
least in the perturbation theory. 

We also note that by the Z^ R-symmetry, the cubic quasi-homogeneous superpotential 
is uniquely singled out. The Yukawa coupling terms ■?/5(l -|- 73)0"'(1 + 73)'?/' with n ^ 1 and 
their conjugates, which may appear in the models other than the A2 model, are not allowed 
by the Z^ R-symmetry. Therefore, we do not need to worry about operator mixings even at 
finite lattice spacing a. 

Unfortunately, however, \D + F\ can be negative. It is easily shown that 7i(-D -|- F)'~fi = 
{D + F)* in the basis 7i = 0-3, 72 = —0-2, implying that every non-real eigenvalue oi D + F 
is doubly degenerated and \D + F\ is real. But the possibility of unpaired real negative 
eigenvalues can not be excluded. Then, |D -|- F| e"*^^ can not be interpreted as a posi- 
tive probability weight. Therefore, our lattice WZ model generally faces the so-called sign 
problem in the standard simulation methods like the hybrid Monte Carlo (HMC) method. 

Sampling configurations — Next we explain how we sample configrations. We utilize the 



Nicolai map 36l-l38| . Using the Nicolai map Eq.(|3]), one can rewrite the expectation value 



of the observable O in the following form, 

(ES^sgn|D + F(0,)|), ' 

where (■■■),, denotes the average over the Gaussian white noise rj. {4>i}i=i,---,N{r]) are the 
solutions of the Nicolai map Eq.(l3]) for a given 77, which contribute to the observable (i.e.|D + 
F(0j)| 7^ 0), and N{ri) denotes the number of the solutions. These solutions exist discretely. 
To see this, suppose a solution (f) + 6(f) that is infinitesimally different from another solution 0, 
\D + F(0)| ^ 0. Then Eq.© means = {Re6(f),ilm6(j)){D + F{(j))) ^ 6(j) = and therefore 
the sigma symbol comes in Eq.Q. 

The expression Eq.(IH]) suggests us a possible interesting simulation method. Suppose 
one succeeds in solving the Nicolai map Eq.(l3]) for a given 77 to obtain the sets {(pi} and 
{sgn|D + F{(j)i)\}. Then, one can simulate the model by observing the numerator and the 
denominator separately, provided that signals for the denominator remain of order 0{1). An 
advantage of this method is that the autocorrelation between samples completely disappears. 
One do not need the updating procedure like the molcular dynamics in the HMC algorithm. 
Of course, one needs to evaluate sgn\D + F(0i)|, but it would not be demading in two 
dimensions for moderate lattice sizes L/a. To solve the Nicolai map numerically, one may 
use the Newton- Raphson algorithm with a globally convergent strategy [39]. A difficulty is 
that one do not know N{r]) a priori. What one may try then is to solve the Nicolai map 
with various initial configurations. 

We sample configurations with the above strategy in this work. For each noise rj, we try 
the Newton-Raphson iterations from 100 initial (p configurations. Since almost all noises 
would be of order 0{1), the solutions should be of order C(l). So we generate these initial 
configurations by the standard normal distribution. And we suppose that the configurations 
obtained in this method exhaust solutions of the Nicolai map. 



We check the quality of these samples by two tests. One is the Witten index 40|]. Note 

that the denominator in Eq.(IH]), A = (X]j=i sgn|D + F(</)j)|)^, is the partition function 

I 
correctly normalized into the Witten index j41j. For instance, in the massive free case 

iy($) = ??i$^/2, the Nicolai map is just a linear equation with a coefficient matrix D + 

m{l — aD/2). Since the determinant of this matrix is positive definite, the solution is 

unique and A indeed reproduces the known result in the continuum theory: (+1)^ = 1. In 



the present cubic case, it must be close to 2 at least in the continuum limit. The other test 
is the Ward identity of the supersjTiimetry Q of the lattice model. For the right hand side of 
Eq.dnj, the on-shell supercharge acts as Qrj = 5S/5ip^, Qtj* = 6S/6ip_. Then {Q{- ■ •)) = 
and the Schwinger-Dyson equations imply the following Ward identities on the lattice: 

{r]{xi) ■ ■■v{xm)v*{yi) ■ ■■v*{yn)) 

m ^ n 

(9) 

J2a'^k=iSx„y^f,^ m = n. 

For instance, the case m, n = 1 reads (Sb) = L^- An ideal case to pass these tests is when 
it happens that X]i=i sgn|D + F\ = 2 over the noises. In this case, the Witten index 2 
is exactly reproduced. In addition, from Eq.([H]), {v ' 'V* ' ' ') reduces to (^ ■ ■ ■ ''7* ■ ■ ■ )r?, and 
the latter reproduces the right hand side of Eq. ([9]) because the noise is the standard normal 
distribution. 

The sampling was carried out at a A = 0.3, on different L x L lattices with L = 
18, 20, 22, 24, 26, 28, 30, 32. We generated 320 noises on each L. The results, shown in 
TABLE [H are classified into following 4 patterns. A: 2 solutions with positive fermion 
determinants. B: 4 solutions, one with negative and the rest with positive determinants. 
C: 1 solution with positive determinant. D: 3 solutions with positive determinants. Since 
almost all the noises belong to A or B where ^jJi sgn|D + F\ = 2, the summation is 
almost 2 over the noises. The Witten index A and the ratio 6 = {{Sb) — L'^)/ Li^ over these 
samples are also shown in TABLE [B From these results, one can estimate the systematic 
error in our simulation method as less than 0.5%. 

Numerical results — With these configrations of the scalar component (f) and the random 
source ?7, one can calculate the various correlation functions in the WZ model. To extract 
the conformal weight of the chiral superfield $, we observe the finite volume scaling of the 
susceptibility of the scalar component 0, 



a^ 



E (<^(^)<^*(o))- (10) 

|x|>3a 



Here we omit the contribution from the correlations at the distances shorter than A~^ = 
a/0.3 ~ 3a, as an improvement of the observable. If the chiral superfield $ really shows the 
scaling behavior with the conformal weight {h, h), and the correlation function {(f){x)(f)*{0)) 



TABLE I: Samples 
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18 


20 22 
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28 


30 


32 
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scales in the IR region, the finite volume scaling of x</. in the continuum limit should read 

/d^x— ^ 

'v F 



X<t> oc 



2h+2/i 



OC V 



l-h-h 



(11) 



where ^ = L^ is the system volume. Using this relation, we read 1 — h — h from the slope 
of the lnx0-lnL^ plots. 

The numerical results are shown in FIG. [1] The error bars show the statistical errors. 
The solid line in FIG. [1] is a fit to the plots by least-square-method and the slope gives 
1 — h — h = 0.660 ±0.011. Our result is consistent with the value expected by the conjecture, 
1 — h — h = 2/3 = 0.666 . . .. This is our main result. 

For a further support of the conjecture, we will also extract the coupling constant K 
of the Gaussian model as follows. In the continuum limit, the WZ model has U{1) R- 
symmetry. Although this chiral symmetry is not broken spontaneously, according to Cole- 
man's theorem 42|, there may appear massive fermions and bosons in the spectrum and 



may decouple in the IR limit, leaving only massless degrees of freedom 43|]. If one writes 
(j) = |0|e*^, the R-symmetry is given by —t- ^ — 2a. Then, the modulus \(f)\ and the new 
fermion x = e^"'^^^'^ip, X = ipe^"'^^^'^ are all singlets and may aquire masses. These degrees of 
freedom may decouple in the IR limit, leaving 6 as low energy degrees of freedom. The IR 
effective action of 6 may be given by 

s,^ = -^ j d?x{vef {e^e + 2Tx), (12) 

where K is an effective coupling constant. Other chiral symmetric terms are irrelevant. 

The value of K will provide us a criterion for the full recovery of A/" = 2 superconformal 
symmetry in the IR limit. The Gaussian model can reproduce M=2 superconformal algebra, 
because both bosonic and fermionic components of the chiral superfield can be constructed 
from the single bosonic field 6*, and the model contains indeed the superconformal stress- 



energy tensor at the two values K = l/127r and S/Air [30|, |3l|- U(l) R-charges of and ifj 
in the WZ model suggest that K = 3/47r is realized in our case. 

To extract the value of the coupling constant K, we again plot the finite volume scaling 
of the susceptibility 

^,^a'J2 {e^'^^'^e-''^'^). (13) 

a::|>3a 

If the above scenario works, in the IR region, the operator e*^ becomes the vertex operator 
with the weight {h^h) = {1/8itK,1/8ttK). Then, one expects xe oc \/0-666..._ ^jj-^g scaling 
dimension is identical with x^- 

The numerical results are shown also in FIG. [1] The solid line is again a fit to the 
plots. The slope is 0.671 ± 0.014 and therefore K = 0.242 ± 0.010, which is consistent with 
K = 3/4tt = 0.238. . .. Thus, our numerical results support the conjecture, including the 
fact c = 1. In addition, since K = 3/47r is the M = 2 supersymmetric point of the Gaussian 
model, this result implies the restoration of all supersymmetries in the IR limit (at long 
distances) without fine tunning. 
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